Orbital order and possible superconductivity in LaNiOs/LaMOs superlattices 
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A hypothetical layered oxide La2NiM06 where Ni02 and MO2 planes alternate along the c-axis 
of ABO3 perovskite lattice is considered theoretically. Here, M denotes a trivalent cation Al, 
Ga,. . . such that MO2 planes are insulating and suppress the c-axis charge transfer. We predict 
that correlated Cg electrons in the Ni02 planes develop a planar x^ — j/^ orbital order driven by the 
reduced dimensionality and further supported by epitaxial strain from the substrate. Low energy 
electronic states can be mapped to a single-band t — t' ~ J model, suggesting favorable conditions 
for high-Tc superconductivity. 
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Despite decades of extensive research, cuprates re- 
main the only compounds to date hosting the high- 
temperature superconductivity (SC). On empirical 
grounds, the key electronic and structural elements that 
support high Tc values are well known ~ no orbital de- 
generacy, spin one- half, quasi two-dimensionality (2D), 
strong antiferromagnetic (AF) correlations. While these 
properties are partially realized in various materials {e.g., 
layered cobaltates) , only cuprates do possess all of them. 

A unique feature of the high- Tc cuprates is the presence 
of an extended doping interval 0.05 ^ (5 < 0.20 where the 
correlated electron maintains its (plane-wave/localized- 
particle) duality, and both fermionic and spin statistics 
may operate in physically relevant energy scales. Mul- 
tifaceted behavior of electrons results in an exotic "nor- 
mal" state of the cuprates with ill-defined quasiparticles, 
pseudogap etc, which condenses into the superconducting 
state below Tc. There are a number of strongly correlated 
metallic oxides [l| based on S = 1/2 3d-ions as Ti'^+ and 
V^^ (both with a single t2g electron), Co^'^ (a t2g hole) 
and Ni'^"'" (closed t2g shell plus one eg electron) that pos- 
sess a low-spin state in octahedral environment. These 
compounds show a great diversity of physical properties 
[l|; however, the mysterious strange-metal phase from 
which anomalous SC may emerge is missing. 

Apart from dimensionality, the orbital degeneracy is 
"to blame" here. Originating from high symmetry of the 
MeOg octahedron - a common building block of both 
pseudocubic and layered perovskites, - the orbital de- 
generacy enlarges the Hilbert space and relaxes kinemat- 
ical constraints on the electron motion. Consequently, 
a fermionic coherency is enhanced and doping induced 
insulator-metal transitions occur without a reference to 
the pseudogap phase. E.g., in Lai_a;Sra;Ti03 the forma- 
tion of a three-band, correlated Fermi-liquid completes 
within just a few percent doping range near x ~ 0.05 [3]. 

The orbital degeneracy strongly reduces AF correla- 
tions (believed to be crucial in cuprate physics), as elec- 
trons are allowed to have parallel spins residing on the 
different orbitals. This leads to competing Ferro- and 



AF-interactions that result in a rich variety of magnetic 
states in S" = 1/2 oxides RTiOa, Na^Co02, Sr2Co04, 
RNiOa, NaNi02, etc. In contrast, spin correlations in 
single-band cuprates are of AF nature exclusively and 
hence strong. 

How to suppress the orbital degeneracy and promote 
cuprate-like physics in other S = 1/2 oxides? In this 
Letter, we suggest and argue theoretically that this goal 
can be achieved in oxide superlattices. Specifically, we 
focus on Ni-based superlattices (see Fig. [T]) which can 
be fabricated usiiig recent advances in oxide heterostruc- 
ture technology ([3|-l5| and references therein). While the 
proposed compound has a pseudocubic ABO3 structure, 
its low-energy electronic states are confined to the Ni02 
planes hence are of quasi-2D nature. A substrate induced 
compression of the NiOg octahedra further stabilizes the 
x^ — y^ orbital. Net effect is a strong enhancement of AF 
correlations among spin one-half electrons residing pre- 
dominantly on a single band of x'^ — y^ symmetry. Thus, 
all the "high- Tc conditions" are perfectly met. More- 



over, the presence of a virtual 3z 



orbital deforms 



the band dispersion compared to cuprates, leading to an 
enhanced next-neighbor hopping t' which is known to 
support higher Tc values [6|. 

The proposed superlattices can be viewed as a layer- 
by-layer "mixture" of a correlated Cg-band metal LaNiOs 
and a band insulator LaMOa. The MO2 planes with a 
trivalent M=A1 or Ga serve here as block layers suppress- 
ing the c-axis hopping. The lattice parameters of LaAlOs 
(~ 3.79A) and LaGaOa (~ 3.89A) are close to that of 
LaNiOs (~ 3.83A); further, Ni02 and MO2 planes have 
the same nominal charges. These factors should result 
in a minimal only structural and electronic mismatch, 
suggesting a stability of La2NiM06 compounds. Yet an- 
other intriguing option is the case of M=Ti (the lattice 
parameter of LaTiOa ~ 3.96A), where Ti02 planes would 
themselves have a spin one-half residing on the t2g or- 
bital. As the t2g shell of Ni'^"'" is full while eg level of 
Ti^^^ is located well above the Fermi energy (roughly at 
lODq ^ 2 eV), the c-axis hoppings will be strongly sup- 
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FIG. 1: (a) Superlattice La2NiM06 with alternating Ni02 
and MO2 planes. MO2 layers suppress the c-axis hopping re- 
sulting in 2D electronic structure. Arrows indicate the c-axis 
compression of the NiOe octahedron imposed by tensile epi- 
taxial strain and supported by Jahn- Teller coupling. (b,c,d) 
Strain induced stretching of the Ni02 planes occurs when 
superlattices with M=A1, Ga, Ti respectively are grown on 
SrTiOa or LaGaOs substrates having large lattice parame- 
ter compared to that of LaNiOa. Expected deformations are 
indicated by arrows. 



pressed again; therefore, Ni02 and Ti02 planes both de- 
velop quasi-2D electronic states of Cg and t2g symmetry, 
respectively. While these states are not mixed by sym- 
metry, the exchange of AF spin fluctuations will lead to 
a sizable interplane coupling, having interesting implica- 
tions for magnetism and possible SC in La2NiTi06. 

The c-axis compression of the NiOe octahedra favor- 
ing x^ — y^ orbital can be imposed by epitaxial strain 
from the substrate with a lattice parameter larger than 
that of LaNiOa, e.g., on LaGaOg or SrTiOs (=i 3.90A). 
The orbital selection by tuning epitaxial strain has been 
demonstrated in Ref. [3| ; it is based on strong Jahn- Teller 
response of the eg orbital on volume conserving distor- 
tions of the oxygen octahedron [3]. In superlattices of 
alternating Ni02 and Ga02 planes (Fig. [T];), strain ef- 
fect should be more efficient as the lattice constants of 
LaGaOa and SrTiOs nearly match. Hence, lattice re- 
laxation is reduced and thicker superlattices La2NiGa06 
can be grown. Further, using the LaGaOa substrate 
(instead of SrTiOa) may help to disentangle the intrin- 
sic physics within the hulk of La2NiM06 superlattice 
from the charge-transfer effects at the interface between 
La2NiM06 and the SrTiOs substrate, induced by differ- 
ent valences of Sr and La [J, |8| . 

Now, we turn to the theoretical examination of our pro- 
posal. The formal valence state in undoped case is Ni'^^ 
in a low-spin configuration ifg^s' ^^"^ ^^^ 4- fold degener- 
acy is specified by spin 5 = 1/2 and orbital pseudospin 
T — 2- Because of strong pd-covalency, wave functions 
are composed of the Ni eg states and a proper combi- 
nation of the oxygen p holes of the same Cg symmetry. 
In other words, the Ni'^^ is shorthand notation for the 



(NiOg/2)'^~ complex with the same - spin doublet ® 
orbital doublet - quantum numbers. Nearest-neighbor 
(NN) hopping matrix within the Ni02 plane, as dictated 
by symmetry, is: 



ta/i — 
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tV3 



TV3 



ct,(3^ {x, z} 



(1) 



Here, {x = x"^ — y^,z = 3z^ — r^} basis is used, =F sign is 
valid for the a and b bond respectively. We assume small 
c-axis hopping rjto ^ to through the block MO2 layers. 

The tensile strain effect is modeled by orbital split- 
ting A = e^ — £2;. The (volume conserving) strain from 
SrTiOa substrate may provide ^ 5% difference between 
the long/short axes of NiOe octahedron, nearly half of 
that in LaMnOs where this leads 0.5 — 1.0 eV orbital 
splitting. Thus, A '^ 0.2 — 0.4 eV might be realistic. 

The 3D nickelate LaNiOa is a correlated two-band 
metal. Nickelates of a smaller-radii rare-earth ions, where 
the bandwidth is reduced due to a stronger GdFeOa-type 
distortion, undergo the insulating state at low tempera- 
ture |9l4l2| and show a peculiar AF order 13[. Whether 
the 2D Ni02 planes of La2NiM06 are insulating or not 
has to be decided by experiment. Physically, the re- 
duction of c-axis hopping and crystal-field splitting both 
should support the insulating state via the orbital dispro- 
portionation phenomenon in correlated systems |14i] . In 
other words, correlations are effectively enhanced when 
the orbital degeneracy is lifted 15[. Based on these ar- 



guments, we consider below the insulating ground state. 
We derived the superexchange Hamiltonian, including 
both intersite (dd) Ni^+-l<ii^+ ->Ni'*+-Ni^+ and charge- 
transfer (CT) [ig Ni3+-02--Ni3+ -^ Ni2+-0-Ni2+ 
processes, along the lines of Ref. llT). The result for a 
bond ij II 7 in the Ni02 plane can be written as 



ijf;) = jo«±.± + if.^I.±)(i ± &^) (h ± &^' 



){^±r^^%±r^^^) 



Pa 



(2) 

with an implied sum over cr = 0, 1 and all combina- 
tions of ±,±. Here, Jo — At^/U , and U stands for 
the Coulomb repulsion on the same orbital. In the c- 
axis, CT part of ([2]) drops out and the rest is multiplied 
by small rj^ . We have used projectors to a singlet and 
triplet state of two Ni^+ S — 1/2 spins: Po,ij — \ — SiSj 



biJ'i 



respectively, and orbital projec- 

.(7) 



e.g., x" — y" and 3z^ — r' 
= c. While the form ^ 



and Pi^ij = 

tors i -I- T^^ and ^ — t^'^' selecting the planar orbital 
in the plane perpendicular to the 7-axis and the direc- 
tional orbital along this axis, 
orbitals, respectively, when 7 
is determined by symmetry, the coefficients K'^'^''-'^ are 
sensitive to the multiplet structure of excited states, as 
usual. Namely, they depend on two parameters Jh/U 
and 2U/{2Apd + Up) characterizing the strength of the 
Hund coupfing and the strength of CT processes [l6|, |l7| , 
resp ectively (explicit expressions will be given elsewhere 
[l8|). In the following, we study the mean-field phase 
diagram depending on these two parameters. 
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FIG. 2: Phase diagram of the spin-orbital Hamiltonian ((2|. 
The c-axis hopping is 0.3to- The labels AAA etc indicate the 
magnetic order, with F (A) denoting parallel (antiparallel) 
orientations of the NN spins in the a, h and c directions, re- 
spectively. Representantive orbital states are sketched. In the 
pure planar-orbital phase (left-bottom), a weak A-coupling 
between the layers is due to the higher order hoppings. Up- 
per right panel; the phase diagram for A = e^ — e^, = 0.5 Jo, 
the gray area indicates the planar-orbital/spin-AAA phase. 
Lower right panel: the phase diagram at 5% hole doping 
(A = 0, J = 0.3to) shows that charge carriers strongly fa- 
vor the planar orbital phase. 



(a) 



zi. 

I 



(b) 



zi. 

I 




12 3 4 5 6 





8=0.2 






1.0 


V 




r 


0.5 


-\^ 




- 


0.0 


- 


/- 


J _ 


0.5 


- V 


/ 


- 




0.5 0.0 0.5 1.0 




0.2-0.1 0.0 0.1 0.2 0.3 0.4 
E/tn 



FIG. 3: The band structure and density of states for the eg 
orbital splitting e^ ~ e^ — 0.25fo at doping levels (a) 100% 
(bare bands), (b) 20%, and (c) 5%. The dotted curves show 



the density of states projected onto the 3z 



orbital. The 



gray areas indicate the occupied states at given S. 



The orbital state |3z^ — r^) cos ^6i + \x^ — y^) sin ^9i 
is specified by the site-dependent orbital angle 9i. We 
assume two sublattices with 6a and 9b in each Ni02 
plane and different spin arrangements and minimize su- 
perexchange (H) in this state. The result is the phase 
diagram in Fig. [2] The antiferromagnetic phase in the 
lower left part consists of pure x'^ — y^ planar orbitals 
(as in cuprates). At higher strength of CT processes, 
this phase changes to the mixed orbital phases AAF, 
AAA, and AFA {9a — 9b evolve gradually from it to 
w 7r/2) with strongly anisotropic spin couplings. Upper 
in-plane ferromagnetic phases have staggered orbital or- 
der 9a = '2tt-9b^ 7r/2. 

Effect of the strain-induced field A on phase diagram 
is dramatic. Even a very modest splitting A makes the 
planar orbital phase a dominant one (Fig. [5]). 

Next, we address the effect of mobile holes (in- 
troduced either by an oxygen-excess or by Sr dop- 
ing). Evaluating ([T]) in our orbital state specified by 
{9i}, we arrive at the effective NN hopping {tap)ij — 
to (2cos6l_ - cos 6*+ =F \/3sin6'+) /4, where 9± = {9i ± 
9j)/2. {E.g., in the planar orbital phase with 9 = tt, 
one finds (i) = 3io/4 as in cuprates.) For a simple esti- 
mate, we add kinetic energy — 2|(<)q| — 2|(i)f,| per hole to 
the superexchange energy. The modified phase diagram 
in Fig. [5] shows that the planar-orbital phase gains the 
largest kinetic energy [13| and quickly spreads with S. 

The spin-exchange constant in the planar orbital phase 
takes the familiar form (at Jh/U = 0) [16J: J = 



(9^2/4) [1/U + 2/{2Apd + Up)]. With t = 3to/4 (hopping 
between the planar orbitals), the prefactor reads 4i^ as 
in cuprates. The value of J is not much sensitive to the 
strength of the Hund coupling, it is about 10% smaller 
at finite Jh/U = 0.15. 

At larger doping, the kinetic energy dominates and the 
spin order is lost. We show in the following that the or- 
bital degeneracy is lifted again and the planar orbital is 
selected by correlated motion of holes. To address this 
regime, we adopt the multiorbital Gutzwiller approxima- 
tion as described, e.g., in Ref. l20l . The hopping matrix 
elements are renormalized tap — >■ Qaptap according to 
the orbital occupations via gap — 2S/y^{2 — na){2 — np) 
(where a,l3 € {x,z}) and orbital splitting is adjusted 
A — > A to minimize energy. The renormalized disper- 
sion Efc is then 



2efc 

to 



{i 



.9+7+ ±Wg_7++ A/to + 12gL7^ (3) 



apart from a fc-independent constant. Here g± = 3gxx i 
g^z, and 7± == ^{cos ka ± cosh). 

Fig. [3] shows bands {sk — M)/io and density of states 
for different dopings at bare orbital splitting A — 0.25to- 
Correlations separate the bands and suppress the or- 
bital mixing, such that the lower band is dominated by 
x^ — y^ orbital and contains only a negligible 3z^ — r^ 
fraction. Correlation-induced increase in orbital split- 
ting (see Fig. 2^) and a reduced mobility of the 3z^ — r^ 
orbital cooperate to lift the orbital degeneracy. The re- 




FIG. 4: (a) The correlation-induced correction to the orbital 
splitting, and (b) the effective t' /t as function of doping 5 for 
different e^ — e^. (c) The mean- field Tc [24I at various levels 
of eg orbital splitting. The dashed-dotted line (e^ — e^ = 5to) 
corresponds to the "cuprate" situation with inactive Zz^ — r^ 
state. 



maining effect of the initially degenerate 3z^ — r^ orbital 



is a deformation of the lower band of : 



y symmetry, 



so that the van Hove singularity stays close to the Fermi 
level. 

The Fermi surface shape can exactly be reproduced 
by an effective t — t' model with the dispersion relation 
— 2t(cos ka + cos kb) — it' cos ka cos kb. The finite t'/t val- 
ues (Fig. \^) are generated solely by the presence of the 
second band, as we started with NN-only hopping. 

Having a single band operating near the Fermi-level, 
we may consider t — t' ~ J model as in cuprates. To see 
effects of strain- induced splitting A (via Sk) on possible 
SC, we estimated the Tc values within the mean-field 
treatment of t — J model (see, e.5., Ref. |211 ). i.e. from 



1 = 2J^ 



ll 



£k -Ml 



■ tanh 



|gfc-Ml 

2Tr 



(4) 



We adopted here J = 0.3t sa 0.2io (a reasonable value 
for the planar orbitals), and e^ is the lower branch of 
Eq. ^ . The resulting T^ curves are shown in Fig. 21 As 
already observed in Fig. [31 the presence of virtual 3z^ — r^ 
orbital deforms the band dispersion and enhances the 
density of states near the Fermi level. The corresponding 
Tc enhancement is clearly visible in Fig. JH 

As a final remark, superlattices with two subsequent 
Ni02 planes, e.g. ../NiNi/Ga/MNi/.. imitating double- 
layer cuprates should be interesting. As the planar or- 
bital ordering is quite robust, a breakdown of quasi 2D 
single-band picture is expected to occur at some larger 
critical number of the subsequent Ni02 planes 23|. 

To conclude, we pointed out that the orbitally- 
nondegenerate spin one-half electronic structure ~ as in 



cuprates - is expected in the Ni-based superlattices. This 
suggests that artificially tailored superlattices may open 
new perspectives for the high-Tc superconductivity. We 
hope that the theoretical expectations for our particu- 
lar proposal - "the double perovskite" La2NiM06 - are 
encouraging enough to motivate experimental efforts. 
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